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ABSTRACT 

Many of the faint submm sources uncovered by deep SCUBA surveys still remain 
unidentified at most, if not all other wavelengths. The most successful method so 
far of obtaining accurate positions and hence allowing more secure identifications 
has been to use centimetre wavelength imaging with the VLA. This has led to a 
tempting approach for obtaining redshift estimates for this population (Carilli & Yun 
1999,2000), which relies on the tight FIR-radio relationship and takes advantage of 
the steep spectral slope in the submm. In this paper we use the submm data from the 
SCUBA Local Universe Galaxy Survey (SLUGS) to estimate the usefulness of, and 
the uncertainties in, the radio-submm redshift estimator. If the submm-radio spectral 
index were correlated with either dust temperature or 850/im luminosity, this method 
could produce biased redshift estimates for 850/zm selected galaxies. We find, however, 
that within SLUGS, these correlations are not significant. The ratio of 850/im/1.4 GHz 
flux was found to decrease with increasing radio and FIR luminosity, and we propose 
that this is due to a component of the dust not associated with recent star formation, 
but which is instead heated to 15-20 K by the general interstellar radiation field. 
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1 INTRODUCTION 

Ever since the discovery of Ultra-luminous Infrared Galax- 
ies (ULIRGs) (Joseph et al. 1985; Sanders et al. 1988) which 
emit up to 99 per cent of their bolometric luminosity in the 
far-IR, there have been suspicions that our optical/UV view 
of the early universe might be biased by large amounts of 
dust obscuration. These were supported by the discovery of 
the Cosmic infrared Background (CIB) (Puget et al. 1996; 
Fixsen et al. 1998), which contains up to twice the luminos- 
ity of the integrated optical/UV background. It thus seems 
likely that much of the star formation activity in the early 
universe takes place in obscured environments. Many deep 
surveys have been conducted with the SCUBA bolometer 
array on the James Clerk Maxwell Telescope (JCMT)* in 
order to uncover the sources responsible for this background. 
Significant numbers of sources, far in excess of no evolution 
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models have been found at both submm (Smail, Ivison & 
Blain 1997; Barger et al. 1998; Hughes et al. 1998; Eales et 
al. 1999) and far-IR (Puget et al. 1999; Kawara et al. 1998) 
wavelengths. 

The next step after finding these sources is identifying 
them and determining their redshifts, thus allowing the true 
star formation history of the universe to be traced. This has 
proved to be a difficult task, due to the large SCUBA beam 
size at 850/im (15 arcsec) and the several arcsec positional 
uncertainty inherent in all of the SCUBA maps (Ivison et 
al. 1998,2000; Downes et al. 2000; Barger, Cowie & Richards 
2000; Eales et al. 2000). The SCUBA sources can sometimes 
be identified with several plausible optical candidates, or in 
other cases they remain optical blank fields, even to the mag- 
nitudes probed by HST and Keck. A possible solution to the 
problem of poor positional accuracy lies in the tight corre- 
lation between thermal FIR emission and synchrotron radio 
emission seen in the local universe (Helou et al. 1985; Con- 
don 1992). The mechanism for this relationship is believed 
to be the massive stars from recent star- formation, which 
both heat the dust to produce the FIR flux and, when they 
explode as supernovae, provide relativistic electrons which 
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constitute the synchrotron radio emission. This means that 
radio observations with the VLA are also sensitive probes 
of star formation, albeit not as sensitive as the submm at 
higher redshifts given the unfavourable K-correction at radio 
wavelengths. For example, Barger, Cowie & Richards (2000) 
have recently targeted 1.4 GHz sources in the deep VLA im- 
ages (Richards 1999) of the Hubble Flanking Fields (HFF) 
with SCUBA and found that the majority of bright (> 6 
mjy) submm sources have radio counterparts. By observing 
the submm sources from the deep surveys at radio frequen- 
cies (particularly 1.4 GHz), more accurate positions can be 
obtained, allowing more secure identifications of the submm 
sources with optical or near IR counterparts. However, even 
with the accurate VLA positions, many SCUBA sources still 
have no detectable optical/IR counterparts, or they are too 
faint for spectroscopy. This creates a great problem in de- 
termining the redshifts for the sources, necessary to trace 
their evolution. 

One possible route for obtaining redshift estimates for 
these objects has been suggested by Carilli & Yun (1999). 
This is based on the FIR-radio correlation plus the break in 
the spectral slope at around 3 mm, where the thermal dust 
emission takes over from the declining synchrotron tail. This 
puts the radio on one side of the spectral break and the FIR 
(or submm in this case) on the other, creating a redshift 
sensitive ratio. Carilli & Yun (1999) used the submm-radio 
spectral index, defined as: 

850 n -r, l ( "->850 \ 

"1.4 = 0.42 x log J 

where Ssso and Si. 4 are the fluxes at 850/im and 1.4 GHz. 

The exact dependence of af'® on redshift is sensitive to 
a few parameters: 

(i) The slope of the thermal Rayleigh- Jeans tail in the 
FIR and submm. This depends on frequency as v 2+l3 where 
the dust emissivity index ((3) is thought to lie between 1 and 
2 (Hildebrand 1983). It is the steepness of this slope which 
potentially makes the indicator sensitive to redshift 

(ii) The slope of the radio synchrotron emission. For star 
forming regions this is assumed to be —0.7 =>- —0.8 (Con- 
don 1992). At higher frequencies this will flatten due to radio 
free-free emission but this is not simple to model. Free-free 
absorption at lower frequencies can also lead to a flatten- 
ing of the radio spectra, the likely effect of this on af 5 f is 
discussed in Carilli & Yun (2000) and they conclude that 
it is unlikely to be a dominant cause of scatter in the rela- 
tionship. In any case the slope at radio frequencies is much 
shallower than in the submm and any uncertainties in it 
have less effect. 

(iii) The temperature of the dust. This determines the 
redshift at which the thermal spectrum turns over. When 
this happens the ratio only depends weakly on redshift, and 
so becomes useless as a redshift indicator. Accurate dust 
temperatures are notoriously difficult to estimate and re- 
quire measurements at many different wavelengths through- 
out the submm and FIR. 

(iv) Contribution by AGN. Galaxies containing a radio- 
loud AGN will have higher radio emission than a purely 
star- forming galaxy and so lower values of a? 5 !*. This could 
lead to ambiguity between high redshift objects with AGN 
and lower redshift starbursts. 

(v) The FIR-radio correlation and its variability with red- 



shift. In using this technique to estimate redshifts, an inher- 
ent assumption is made that the FIR-radio correlation does 
not vary with redshift. Possible causes of such a change could 
be a variation in magnetic field strengths in galaxies in the 
past, or changes in the dust mass opacity coefficient which 
relates the FIR/submm emission to the mass of dust. The 
latter would depend on the dust composition and grain sizes 
so there is scope for it to have been different at the epoch of 
galaxy formation, when metallicities would have been lower. 

(vi) Linked to (v) is the possibility that the FIR-radio, or 
more importantly, the submm-radio relationship itself de- 
pends on some other galaxy property such as dust temper- 
ature, luminosity etc. 

The possible impact of points i-vi above will be discussed 
later. 

Several authors have already used this technique to de- 
termine redshift estimates or limits for sources discovered in 
deep submm surveys with SCUBA (Hughes et al. 1998; Lilly 
et al. 1999; Smail et al. 1999; Barger et al. 2000; Eales et 
al. 2000). This paper uses data from the SCUBA Local Uni- 
verse Galaxy Survey (SLUGS) to examine the usefulness of 
the submm-radio ratio as a redshift indicator, to try to nor- 
malise it with respect to the local universe and to estimate 
the uncertainties in it. 

We will use H = 75 km s _1 Mpc -1 and q = 0.5 
throughout. 



2 CONSTRAINING THE REDSHIFT 
INDICATOR 

The low redshift data is taken from the SLUGS survey 
(Dunne et al. 2000) which is a complete sample of 104 
galaxies selected from the IRAS Bright Galaxy Sample. 
The sample includes mainly infra-red bright objects (La T = 
7x 10 9 -2x 10 12 L ra ), some ULIRGS as well as many 'more 
normal' galaxies. The objects were mapped at 850/im with 
SCUBA and so should provide reliable submillimetre fluxes. 
We fitted a single-temperature dust spectral energy distribu- 
tion (SED) to the measured fluxes for each galaxy, allowing 
the dust emissivity index (/3) as well as the temperature to 
vary to find the best fit. While the assumption of a single 
temperature is probably wrong, and indeed the low values 
of (5 we derive suggest that there is a cold dust component, 
the SEDs are still a good empirical fit to the data, which is 
all that we require for the present analysis. 

The original Carilli- Yun models for the behaviour of the 
radio-submm spectral index with redshift were simple two 
power-law models, along with some empirical SED data for 
2 local galaxies. They have since modified their original esti- 
mator by using data from nearby galaxies (17 objects) to try 
to provide a normalisation and to estimate the scatter in the 
submm-radio ratio with redshift (Carilli & Yun 2000). We 
have independently undertaken a similar analysis but using 
the full SLUGS data set (104 objects), thereby providing a 
more statistically valid measure of the relationship and its 
scatter. 

We have taken the real spectral energy distributions of 
the 104 SLUGS galaxies as determined from the IRAS 60 
and 100/im fluxes and the SCUBA 850/im data, and the 
1.4 GHz radio flux (Condon et al. 1990). The measured ra- 
dio spectral indices were used where available (for about 
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Table 1. la uncertainties on Z mB< i 



^med 


±<T Z 


1.0 


0.28 


2.0 


0.44 


3.0 


0.63 


4.0 


0.94 


5.0 


1.25 


6.0 


+3.13 




-1.44 



Uncertainties based on scatter in low redshift data only, 
there will be additional uncertainty when the flux errors 
of the high-z sources are accounted for. 

half of the sample) and the median a ra dio was found to be 
—0.7. This was used in the SEDs of those galaxies with no 
measured a ra< n - If the change in each of these SEDs with 
redshift is plotted, we get 104 separate curves from which 
the median can be selected. The uncertainty in this method 
(at any given redshift) can be estimated from the spread in 
the 104 curves. In Figure 1, the thick solid line shows the 
median curve with the upper and lower solid curves being 
those ranked 16 per cent from the top and bottom respec- 
tively, to give the ±1<7 errors (defined as the range of a? 5 ! 1 in 
which 68 per cent of the lines lie). So at any given redshift, 
the horizontal distance between the thin lines provides an 
estimate of the ±lcr uncertainty on the redshift provided by 
the median curve. The rms scatter in af 2 is almost constant 
with redshift, decreasing very slightly from 0.097 at z = to 
0.080 at z = 6. A rough guide to the redshift uncertainties 
at various redshifts are given in Table 1. The median curve 
can be approximated by the expression 

z = 0.551 - 6.652a + 25.57a 2 - 30.56a 3 + 13.75a 4 

Our estimator starts to loose its effectiveness at redshifts 
greater than about 4-5, where it starts to turn over due 
to the peak in the thermal dust spectrum. This can be seen 
from the increase in a z with redshift (Table 1). Also, as men- 
tioned in Carilli & Yun (1999) at z > 6 the radio emission 
begins to become quenched by inverse Compton losses off 
the microwave background. The dashed line in Fig. 1 shows 
the new Carilli- Yun estimator (with ±la uncertainties) us- 
ing the data from 17 local galaxies (Carilli & Yun 2000). It 
can be seen that their mean estimator lies outside the la- 
uncertainties from our analysis, however the upper la curve 
from theirs is consistent with our lower la limit, with agree- 
ment being better at lower redshifts (< 2). The differences 
in redshift estimates from the two indicators are listed in 
Table 2 for various values of af 2- 

There are two main differences between the present re- 
sults and those in Carilli & Yun (2000). The first is in the 
zero redshift normalisation in the two samples. The mean 
value of a? 5 4° at z = in Carilli & Yun (2000) is w 
with an rms scatter of 0.14, while for our sample the mean 
af 5 ° = 0.18 at z = 0, with a lower rms scatter of 0.097. 
The shape of the Carilli- Yun estimator is also different due 
to the SED shapes of the galaxies in their sample (which 
in general had higher values of /3 than ours). These differ- 
ences are partly due to our much larger sample size, but 
there are also discrepancies in the submillimetre fluxes of 
the individual galaxies which are present in both the sam- 
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Figure 1. The median (thick solid) and ± la (thin solid) esti- 
mators using the SLUGS low redshift data. The recent Carilli & 
Yun (2000) estimator is shown (dashed) with ±1<t uncertainties 
for comparison. The differences in redshift estimates between this 
work and the Carilli- Yun estimator at various values of a^ 5 ® are 
given in Table 2. 

pies used by ourselves and Carilli & Yun (2000). Most of 
their submm data comes from a sample of 19 bright IRAS 
galaxies (Lisenfeld, Isaak & Hills 2000) . Seven of the galax- 
ies in this sample were also observed by us, and in several 
cases the fluxes we derived were as much as a factor of two 
higher. This was attributed to Lisenfeld et al. missing flux 
from extended objects (Lisenfeld, private communication), 
and has been corrected in the revised version of the Lisen- 
feld et al. paper. The consequences of using under-estimated 
submm fluxes would be to i) lower the mean af 5 ° and ii) in- 
crease the value of /3 determined in single temperature fits 
to the 60, 100 and 850/im fluxes, which explains the different 
shape of the Carilli- Yun estimator. We also note a possible 
error in the 1.4 GHz fluxes used for Zw049.057 and NGC 
5936 by Carilli & Yun (2000). An erratum has since been 
published to accompany Carilli & Yun (2000) and this is 
discussed further in Section 5. 

If we now place the handful of SCUBA sources from 
the deep surveys which have both spectroscopic redshifts 
and radio fluxes (Smail et al. 1999 and references therein; 
Eales et al. 1999; Lilly et al. 1999; Barger et al. 2000; Eales 
et al. 2000; Ivison et al. 2000) on the redshift estimator we 
have created from our observed SEDs (Fig. 2), we see that 
the agreement is satisfactory given the uncertainties. The 
object at z = 2.8 is SMM02399-0136 from Ivison et al. (1998) 
and is known to harbour an AGN, which explains why it 
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Figure 2. The median (thick solid) and ± lcr (thin solid) estima- 
tors using the SLUGS low redshift data. Data points are sources 
from the deep submm surveys which have secure radio fluxes and 
spectroscopic redshifts. Triangles are from the Canada-UK deep 
survey (Eales ct al. 1999; Lilly et al. 1999; Eales et al. 2000). 
The point at z = 1.15 is from Eales et al. 2000 and has fluxes 
of 5850 = 3.3mjy, SsGHz = 53.6/iJy with a ra ^ io = —0.3. Stars 
are from the Cluster survey (Ivison et al. 1998; Smail ct al. 2000; 
Ivison et al. 2000). This includes a z = 2.22 galaxy which has 
a radio upper limit. Circles are from the SCUBA observations 
of radio sources in the HFF (Barger et al. 2000). The point at 
z = 2.8, which lies well below all the curves, is SMM02399-0136 
from Ivison ct al. (1998) and is believed to be an AGN. 



Table 2. A comparison of estimated median redshifts and ranges 
of redshift for the relationship in this work and that of Carill & 
Yun (2000) 



850 
"1.4 


SLUGS 


^rned 

CY 


Z+" - z- a 
SLUGS 


z +a _ z -a 

CY 


0.3 
0.5 
0.7 
0.9 
1.1 


0.19 
0.64 
1.20 
2.12 


0.49 
0.96 
1.65 
2.95 
6.04 


0.07-0.34 
0.42-0.82 
0.93-1.54 
1.73-2.62 
2.92-4.46 


0.23-0.81 
0.60-1.42 
1.12-2.45 
1.93-4.44 
3.63-> 7 



Column(l)- Value of spectral index; Co/umn^-redshift 
estimate using median line from SLUGS; Column(3)- 
as for (2) but using Carilli & Yun (2000) estimator; 
column(4)-range of redshifts between ± la curves from 
SLUGS; Column(5)-&s for (4) but using CY (2000). 



lies below the predicted lines. Given that sources may have 
AGN activity but cannot be confirmed as such, this redshift 
estimator can give only reliable lower limits. Additionally, if 
taken with another redshift estimation, such as photometric 
redshift, it could help to identify the presence of AGN (a 
point first noted by Carilli & Yun 1999). 



3 THE FIR AND SUBMM VS. RADIO 
PROPERTIES 

The SLUGS sample provides an ideal test for the varia- 
tion of the FIR-radio and submm-radio relationships with 
galaxy properties, such as dust temperatures, luminosities 
etc. which may affect the reliability of af 5 " as a redshift 
indicator, and which is also interesting in its own right. 
Figure 3 shows the FIR-radio correlation for the SLUGS 
galaxies along with the 850pm-1.4 GHz relationship. Both 
are highly significant (Table 3 lists all the correlation coef- 
ficients in this discussion) but the scatter in the FIR-radio 
correlation is smaller (0.06 (15%)) than that which exists 
for the submm-radio (0.097 (25%)). This suggests that it 
is actually the FIR which is physically related to the radio 
rather than the submm. There is also a difference in slope 
between the two, which we will return to later. The FIR- 
radio correlation in the local universe does not depend on 
galaxy luminosity (except at very low Lfi r ) or dust temper- 
ature, something which is confirmed by this sample. How- 
ever if the submm-radio relationship is tested in this way by 
plotting af 5 4 against Lfi r , L 14 and L 850 (Figure 4) then it 
is clear that there is a strong dependence of af 5 ! 1 on both 
Lfir and L1.4 but not on Lsso. The link with Lfi r was noted 
by Carilli & Yun (2000) but the trend is even more evident 
with radio luminosity. The most likely explanation for these 
correlations is that both Lfi r and L1.4 are very sensitive to 
recent star formation in a galaxy but Lsso is not as sensitive, 
as a significant fraction of the 850pm flux may be produced 
by colder dust heated by the general interstellar radiation 
field (ISRF). Therefore Lsso will not change by as great a 
fraction as Lfi r with increasing star formation. This explains 
the larger scatter and steeper slope in the 850pm-1.4 GHz 
correlation. 

One of the two crucial relations for determining whether 
af 5 4 is a biased redshift indicator is that between a? 5 ] 1 and 
Lsso- Since SLUGS was selected at 60pm, it is quite possi- 
ble that the 850pm luminosities of SCUBA-selected galaxies 
are distinctly different from those in SLUGS. If there were 
also a significant relationship between af 5 4 and Lsso, the 
calibration of the Carilli- Yun method using SLUGS galax- 
ies would lead to biased redshift estimates. Fig. 4(c) shows, 
however, that there is no significant correlation. To exam- 
ine this, the dot-dash line on Fig. 4(c) shows the lower end 
of the rest-frame 850pm luminosity range for a galaxy with 
5850 = 4mJy at redshifts of > 2 (f3 = 2 and T d = 50 K, 
lower (3 and Td give slightly higher rest-frame Lsso for the 
same flux and redshift). If we look at the values of af 5 ! 1 for 
the SLUGS galaxies, which correspond to the luminosities 
of the high redshift objects (i.e. the crosses to the right of 
the dot-dash line in Fig. 4(c)), they do not display a lower 
mean value than that for the whole sample. Given how little 
we know about the properties of the high-z SCUBA galax- 
ies, and whether they are related to galaxies of the same 
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(a) FIR-radio correlation 



(b) 850/im-radio correlation 



Figure 3. The FIR— radio relation and the 850/im-radio relation for the SLUGS galaxies. The FIR-radio correlation is the tighter one 
implying that this is where the physical connection lies. 



luminosity at low redshift (density evolution) or to galaxies 
of lower luminosity (luminosity evolution), we do not feel 
that there is a need for a luminosity dependent indicator. 
This is further supported by Fig. 2, showing a relatively 
good agreement of the indicator with the high- 2 data, and 
in particular, there is no systematic bias from the line in one 
direction, which may be indicative of changes in af 5 !* with 
redshift /luminosity. 



The second important relationship is that between a'^ 
and dust temperature, since an 850/im selected sample may 
well have a different mean dust temperature than one se- 
lected at 60/im. It has also been suggested by Blain (1999) 
that there should be a dependence of af 5 2 on dust temper- 
ature, and that this produces a degeneracy in the redshift 
indicator whereby a hot galaxy at high redshift may be in- 
distinguishable from a colder galaxy at lower redshift. We 
have tested for a temperature dependence by plotting af 5 ^ 
against dust temperature, and we find no significant cor- 
relation (Figure 5(a)). The slight correlation seen between 
a? 5 !" and dust emissivity index [3 (Figure 5(b)) is probably 
linked to the relationship of af 5 f with L1.4 and Lr t . Despite 
the significant correlation given by the statistic to this re- 
lationship, it appears to only hold for the lower values of j3 
(0.9-1.4). This could be because a low /3 (when produced 
by fits to 60,100 and 850/im points only) is another possible 
indicator of the fraction of 850/im emission produced by the 
ISRF, rather than directly by star forming regions. Evidence 
for this comes from observations of a sample of optically se- 
lected galaxies presently under way with SCUBA and also 
from studies of NGC 891 and the Milky Way (Alton et al. 
1998; Masi et al. 1995), where the SEDs using only 60, 100 
and 850/im fluxes give a very low f3 of ~ 0.7. These galaxies 
are known to have large fractions of cold dust at T < 20 K, 
however this does not lead to a dependence of a? 2 on the 
fitted dust temperature because that is determined by the 
FIR fluxes which are dominated by the warmer dust. 



Table 3. Correlation coefficients. 



V 
L\A 


X 
Ltir 


0.92 


significance 
1.5 X lO" 20 


S.D 

9a 


L1.4 


i-850 


0.89 


1.1 X lO" 19 


9a 


850 
tt 1.4 


£fir 


-0.43 


1.3 x 10~ 5 


4a 


850 
tt 1.4 


L1.4 


-0.61 


5.4 x 10" 10 


6a 


850 
tt 1.4 


i<850 


-0.22 


0.03 


2a 


850 
tt 1.4 


Tdust 


-0.22 


0.03 


2a 


850 
a 1.4 


13 


-0.42 


2 x 10" 6 


4a 



notes- Co/«mn(3/'-Spearnian rank correlation coefficient; 
Co/jimn^-probability that x and y are unrelated (i.e. of 
accepting the null hypothesis); Cokimn(5,)-standard devi- 
ation at which the null hypothesis is rejected. 

4 DISCUSSION 

The submm-radio spectral index appears to provide a mod- 
erately satisfactory redshift indicator out to redshifts of 
~ 4 — 5 where it ceases to be as sensitive to redshift. The 
uncertainties in the estimated redshifts are obtainable using 
the spread of the local SEDs, and range from < ±0.3 at 
z < 1 to > ±1 at z > 4. Sources which harbour radio loud 
AGN will not follow this pattern as they will lie below the 
indicators, leading to AGN sources which are not recognised 
as such being given misleadingly low redshifts. 

Let us now return to the assumptions and limitations 
discussed in section 1, and apply them to our data sets. 

(i) The submillimetre emissivity index has a mean value 
of 1.3 for the SLUGS sample, as determined by single tem- 
perature fits to the 60, 100 and 850/im fluxes. As discussed 
at length in Dunne et al. (2000) , if a colder dust component 
is present this may not represent the true value of /?, the 
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(a) Dependence of af 5 ! 1 on Lfi r , first noted by Car- 
illi & Yun (2000).' 
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(b) A stronger correlation is seen with L1.4 
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(c) af 5 4 shows no significant correlation with Lgso 



Figure 4. Relationships between a^'® and luminosities for the SLUGS galaxies. The strongest connection is with radio luminosity while 
there is no strong dependence on L850- The dot-dash line in (c) shows a rough lower limit to the rest-frame 850/im luminosity of a 
deep SCUBA source with Ssso = 4mjy, z > 2, /3 = 2 and T<j = 50 K. Lower f3 or T<j will produce slightly higher luminosities. The flat 
flux density-redshift behaviour at 850/im means that sources brighter than 4mjy will have Lgso higher than the line, and vice-versa for 
fainter ones. 



Table 4. Parameters used in 'grey' cold component models. 



Model 


T w (K) 


Tc (K) 


N c /N w 


P 


1 (thick dashed) 


40 


20 


5 


2 


2 (thin dashed) 


40 


20 


10 


2 


3 (thin dotted) 


40 


20 


40 


2 


4 (thick dotted) 


30 


15 


10 


2 



Column(l)- Refers to model in Fig 6; Column(2)- 
Temperature of warm dust component; Column(3)- 
Temperature of cold dust component; Column(4)- Ratio 
of cold dust mass to warm; Column(5)- Value of emissiv- 
ity index assumed 

true value being higher than that determined using a single 
component model. If the true value of /3 were nearer to 2 this 
would increase the rate of change of the ratio with redshift 



as the submm flux in the Rayleigh- Jeans tail oc v 2+l3 (see 
(hi)). 

(ii) The radio spectral index is also somewhat uncertain, 
but being flatter it has much less impact than the submm 
slope. However, if high-redshift dust sources have systemat- 
ically different spectral indices than those at low redshift, 
this would lead to a bias. The range of measured indices for 
the SLUGS sample is ~ —0.2 to —1.0. The situation at high 
redshifts is harder to determine as the value of a ra( uo for 
[i3y sources may depend on the selection frequency, with 
sources selected at v > 5 GHz likely to have flatter radio 
spectra than those selected at 1.4 GHz (Richards 1999). For 
example, the sources from the Canada-UK SCUBA survey 
have 5 GHz radio fluxes from Fomalont (1991) and an aver- 
age a ra dio = —0.38, while SCUBA observations of 1.4 GHz 
sources in the Hubble Flanking Fields were associated with 
radio sources with steeper spectra, a ra d ~ —0.7 (Barger et 
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dust temperature (K) 



(a) No significant correlation is found between a* 5 ^ 
and the fitted dust temperature 



+ t 



emissivity index (fl) 



(b) Possible slight dependency of of 5 !* on emissiv- 
ity index /3, see text for discussion 



Figure 5. Variation of af 5 ^ with fitted dust temperature and emissivity index 
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Redshift 

Figure 6. Median and lcr estimators (solid lines) along with 
models for two-component dust SEDs. Parameters are detailed 
in Table 4: T w = 40 K, T c = 20 K, N c /N w = 5 (thick dashed); 
T w = 40 K, T c = 20 K, N c /N w = 10 (thin dashed); T w = 40 K, 
T c = 20 K, N c /N w = 40 (thin dash-dotted); T w = 30 K, T c = 15 
K, N c /N w = 10 (thick dash-dotted). 

al. 2000). This is most probably only a selection effect, as 
only relatively flat spectrum sources would be detected at 5 
GHz. In general, it is preferable to survey at 1.4 GHz as this 
gives the best sensitivity to high redshift objects. 

(iii) The mean dust temperature in the local galaxies 
could be over-estimated. If a colder component is present 



then it will produce a flattening of the thermal dust spec- 
trum at longer wavelengths when compared to warm dust, 
thus lowering the redshift at which the estimator loses its 
sensitivity to redshift. We can model the implications of 
points (i) and (iii) as they are inter-related. If there is a 
cold dust component in the local galaxies, the true value of 
(3 is likely to be nearer 2 and the mean dust temperature 
lower. If we make assumptions about the temperatures of 
the two components and the relative masses in each (and 
that /3 = 2), we can then produce a new indicator. Fig- 
ure 6 shows these 'grey' indicators for a few values of the 
various parameters which reflect current observational evi- 
dence (Frayer et al. 1999; Alton et al. 1998; Calzetti et al. 
2000). The parameters used are given in Table 4, and all 
'grey' models were normalised at zero redshift to have the 
mean value for the sample (0.18). There is generally no great 
difference until z ~ 2, and at this point the different param- 
eters assumed start to make a larger impact. In general, we 
might expect that the galaxies detected by both SCUBA 
and the VLA would have T w > 30 K, since it is the recent 
star formation which the VLA is sensitive to (and which also 
produces higher dust temperatures). This makes the upper 
(T w — 40 K) curves more plausible although their shape is 
still sensitive to the mass fraction and temperature of the 
cold dust (colder dust or more cold dust will cause the in- 
dicator to turn over earlier and flatten more), and is also 
dependent on the value of (5. If in reality, /3 lies somewhere 
between 1.3 — 2, then this will make the 'grey' models more 
like the median one at low redshifts. Figure 6 suggests that 
the uncertainties in whether there is actually a single dust 
temperature, or in the true value of /3, do not add much to 
the uncertainty we have derived from the single-temperature 
fits. Currently, it is not possible to be more specific about 
the nature of any cold dust components as there are only a 
handful of local galaxies with enough FIR and submm fluxes 
to make a decomposition of the SED feasible. Our knowl- 
edge of dust properties in local galaxies should improve in 
the near future (Dunne et al. 2000). 

(iv) Contributions to the radio flux by AGN will produce 
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misleading redshift estimates if the object is not recognised 
as such. Since the fraction of deep SCUBA sources harbour- 
ing radio-loud AGN is still not well determined, the estima- 
tor must be treated as a statistical tool rather than a redshift 
indicator for individual objects. The redshifts given by the 
upper la curve should be treated as a robust lower limit. 

(v) A variation in the FIR-radio correlation with redshift 
is quite possible if either magnetic fields or dust properties 
were different in the past, although currently this is difficult 
to test. Condon (1992) does however point out that the FIR- 
radio relation in the local universe holds over a large range 
of magnetic field strengths. The tendency of the galaxies 
observed so far in both the radio and submillimetre at high 
redshifts, to lie in agreement with the estimator (Fig. 2) 
suggests that there has been no dramatic evolution in the 
underlying FIR-radio relation, although more data is needed 
to fully investigate this. 



5 CONCLUSIONS 

We have used the submm data from a large, complete sam- 
ple of local IRAS galaxies, along with complementary ra- 
dio data, to define a median redshift estimator using the 
change in spectral slope between the submm and the radio 
at ~ 3mm. The scatter in the data has been used to provide 
la uncertainties on the relationship. The limited number of 
submm sources from the deep surveys with both radio fluxes 
and spectroscopic redshifts generally agree with the estima- 
tor (within the uncertainties), except for one object which 
is known to be an AGN (Fig 2) . The estimator is useful in a 
statistical sense rather than for predicting the redshifts of in- 
dividual objects as there are many uncertainties, especially 
the possible contribution to the radio flux by AGN which 
would lower the data relative to the models thus leading 
to an under-estimated redshift. Our redshift indicator dif- 
fers from the recent one of Carilli & Yun (2000) in terms of 
shape and normalisation. The difference is significant, par- 
ticularly at redshifts > 2 (Fig 1). The discrepancy can be 
attributed primarily to the larger sample size used in this 
work and also, in part, to the under-estimation of submm 
fluxes used in CY 2000. An erratum to CY 2000 has since 
been published which produces a revised redshift estimator 
using the corrected Lisenfeld et al. (2000) submm fluxes. 
This removes ~ 30 per cent of the discrepancy between the 
estimators, and we postulate that the remaining difference 
is a result of sample selection (the Lisenfeld sample is not a 
complete sample and contains a higher proportion of radio 
bright objects compared to SLUGS), and the smaller sample 
size (17 objects compared to 104 in SLUGS). 

The submm-radio spectral index decreases as radio and 
FIR luminosity increase, but shows no strong trend with 
850/im luminosity (Fig. 4). The correlations are believed to 
be due to the sensitivity of the non-thermal radio and ther- 
mal FIR emission to recent star formation, while a signifi- 
cant fraction of the 850/im flux may be due to colder dust, 
and not linked to the star formation rate. The absence of a 
correlation between af 5 2 and 850/im luminosity means that 
any differences in 1/850 between SLUGS and the SCUBA- 
selected galaxies will not lead to biased redshift estimates. 
Importantly, we also find no significant evidence for a sys- 
tematic variation of af 5 f with dust temperature (Fig. 5(a)) 



implying that the redshift-temperature degeneracy hypoth- 
esised by Blain (1999) does not play a dominant part in the 
scatter of the af 5 f - redshift relation. 

We have investigated the impact of neglecting a possi- 
ble colder dust component in the SEDs of the local galaxies 
(Fig 6). While sensitive to the parameters assumed, in gen- 
eral the differences are not very significant for most plausible 
models and the agreement is still within the la uncertainties 
at lower redshifts (< 2). If cold components are present but 
unaccounted for, using the median line in Fig 1 to estimate 
redshifts will most likely lead to an over-estimate of z. 

We would like to thank Chris Carilli and Ute Lisenfeld 
for useful discussions. The support of PPARC is gratefully 
aknowledged by L. Dunne, S. Eales and D. Clements. 
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